Quasi-particle Statistics and Braiding from Ground State Entanglement 
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Topologically ordered phases are gapped states, defined by the properties of excitations when 
taken around one another. Here we demonstrate a method to extract the statistics and braiding 
of excitations, given just the set of ground-state wave functions on a torus. This is achieved by 
studying the Topological Entanglement Entropy (TEE) on partitioning the torus into two cylinders. 
In this setting, general considerations dictate that the TEE generally differs from that in trivial 
partitions and depends on the chosen ground state. Central to our scheme is the identification of 
ground states with minimum entanglement entropy, which refiect the quasi-particle excitations of 
the topological phase. The transformation of these states allows for a determination of the modular 
S and U matrices which encode quasi-particle properties. We demonstrate our method by extracting 
the modular S matrix of a SU(2) spin symmetric chiral spin liquid phase using a Monte Carlo scheme 
to calculate TEE, and prove that the quasi-particles obey semionic statistics. This method offers a 
route to a nearly complete determination of the topological order in certain cases. 
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I. INTRODUCTION 

Topologically ordered states are exotic gapped quan- 
tum phases of matter that lie beyond the Landau sym- 
metry breaking paradignii. Well known examples include 
fractional quantum Hall states, gapped quantum spin liq- 
uids and quantum dimer models^"— . These phases are not 
characterized by correlations or local order parameters 
but rather by long range entanglement in their ground- 
state wave functions^. Emergent excitations in these 
phases can carry fractional quantum numbers and, in two 
dimensions, realize non-trivial statistics on exchanging 
identical particles or on taking one particle type around 
another—. Two important hallmarks of topologically or- 
dered phases are ground-state degeneracy in a space with 
non-zero genus and a finite topological entanglement en- 
tropy (TEE)ii"— . In this paper we show that combined 
together, these two properties can be used to extract key 
aspects of a topological phase such as braiding rules and 
topological spin of quasi-particlea^^. 

Taking a more general perspective, we note that a 
gapped phase of matter remains stable as long as the 
excitation gap between the many body ground state(s) 
and the first excited state remains non-zero. This im- 
plies that all 'universal properties' associated with the 
gapped phase, which we define as the properties that re- 
main invariant with respect to the change of the underly- 
ing Hamiltonian while maintaining the gap, must be en- 
coded in the ground state wavefunction(s). The question 
of extracting such universal properties solely from the 
ground state wave-function(s) becomes especially inter- 
esting for two-dimensional topologically ordered phases. 
In this case, the universal properties are related to the 
anyonic character of gapped excitations, such as braid- 
ing rules of elementary quasi-particles. The preceding 
argument seems to employ that in a topological ordered 
phase, all robust properties of gapped excitations are en- 



coded in the ground state (s) itself. Such a point of view 
was taken by X.-G. Wen in Refi^ where the notion of 
non-abelian berry phase was introduced to extract the 
braiding and statistics of quasi-particles. However, the 
idea in Refi^ requires one to have access to an infinite set 
of ground-states labeled by a continuous parameter, and 
is difficult to implement. 

Recently, it was found that the ground-state entangle- 
ment entropy of a two dimensional topologically ordered 
phase in a disk-shaped region A with a smooth bound- 
ary of length L takes the form Sa = olL — 7, where 
the universal constant 7 is the TEEi^ "— . W hen, The 
constant 7 equals log(-D) where D = y^rff is the "to- 
tal quantum dimension" associated with the topological 
phase while di is the quantum dimension of I'th quasipar- 
ticle type. Unfortunately, the total quantum dimension 
D only provides a partial characterization of topological 
order (for example, two distinct topological phases can 
have same value oi D). A natural question arises whether 
the quantum entanglement can be used to extract indi- 
vidual quantum dimensions and perhaps, the anyonic 
braiding and statistics associated with the quasiparticles? 
As we show in this paper, the answer to this question is 
yes. We will only require the degenerate set of ground 
state wavefunctions on the torus. 

Recall that the ground state has topological degener- 
acy when the system is defined on a topologically non- 
trivial manifold. It is generally assumed that TEE is a 
quantity solely determined by the total quantum dimen- 
sion D of the underlying topological theory as 7 = log D 
independent of which ground state it is being calculated 
for. However, this holds true only when the boundary of 
the region A consists of topologically trivial closed loops. 
If the boundary of region A is non-contractible, for exam- 
ple if one divides the torus into a pair of cylinders, gener- 
ically the entanglement entropy is different for different 
ground states. Indeed as shown in Ref. [l^ for a class of 
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topological field theories, the TEE depends on the par- 
ticular linear combination of the ground states when the 
boundary of region A contains non-contractible loops. 

In this paper, we first elaborate on the ground-state de- 
pendence of entanglement entropy, focusing on the case 
of a partition of torus into two cylinders. We present an 
argument based on the strong subadditivity property of 
quantum information and show that the TEE per con- 
nected boundary is not identical to that for a trivial bi- 
partition, such as a disc cut out of the torus, where TEE 
is logZ^. This is illustrated as an 'uncertainty' relation, 
between entropies for two different cylindrical biparti- 
tions of the torus. 

We also demonstrate this ground-state dependence nu- 
merically, and calculate the entanglement entropy for the 
chiral spin liquidi^ (CSL) wave function as different lin- 
ear superpositions of the two ground states (Figure [3]) 
with the Variational Monte Carlo (VMC) methodil^i^. 
The physical origin of the ground-state dependence of 
TEE is made explicit by studying a Z2 toric code model^. 
We introduce the notion of minimum entropy states 
(MESs), namely the ground states with minimal entan- 
glement entropy (or maximal TEE, since the TEE always 
reduces the entropy) for a given bipartition. These states 
can be identified with the quasi-particles of the topologi- 
cal phase and generated by insertion of the quasi-particles 
into the cycle enclosed by region A. For a generic lattice 
wave function with finite correlation length, such as the 
CSL wave functions we study, a nonlocal measurement 
like TEE is essentially to identify this basis of MESs. 

Having established the dependence of the TEE on the 
ground states, we detail a procedure that uses this de- 
pendence to extract the key properties of quasi-particle 
excitations by determining modular S and U matrices - 
a vital characteristic of topological orde r^i^^i^"" — . Unlike 
the TEE, which failed to distinguish topological phases 
with the same total quantum dimension, modular S and 
U matrices provide a finer characterization of topological 
order. In fact, they fully determine the quasi-particle self 
and mutual statistics as well as the individual quantum 
dimensions of particles in a topological phase. In partic- 
ular, the element Sij of the modular S matrix determines 
the mutual statistics of i'th quasiparticle with respect to 
the j'th quasiparticle while the element Uu of (diagonal) 
U matrix determines the self- statistics ('topological spin') 
of the i'th quasiparticle. This procedure requires just the 
set of ground states, although the results pertain to the 
braiding and fusing of gapped excitations. The basic idea 
is to relate MESs for different entanglement bipartitions 
of the torus. The MESs, which reflect quasi-particle ex- 
citations, are determined using TEE. As an application, 
we extract the modular S matrix of an SU{2) symmetric 
CSL, a lattice equivalence to a = 1/2 Laughlin state, 
through TEE calculated with a recently developed VMC 
scheme. For illustrative purposes, we discuss in the Ap- 
pendix how our algorithm applies to Kitaev's toric code 
model^, a zero correlation length phase, and extract the 
modular S matrix and U matrices. 



At a practical level, the procedure outlined here sug- 
gests that entanglement entropy could be used to numer- 
ically diagnose details of topological order beyond the 
total quantum dimensio n^^i^'^" — , which is a single num- 
ber susceptible to numerical error. An elegant differ- 
ent approach to a more complete identiflcation of topo- 
logical order is through the study of the entanglement 
spectruni^l. However we note that requires the existence 
of edge states and may not be applicable for topological 
phases like the Z2 spin liquid. Furthermore, it is possi- 
ble to compute TEE using Monte Carlo techniques on 
relatively larger system o-^^i^^ , as also done in this paper, 
where the entanglement spectrum is not currently avail- 
able. 



II. GROUND STATE DEPENDENCE OF 
TOPOLOGICAL ENTANGLEMENT ENTROPY 

Given a normalized wave function |$) and a partition 
of the system into subsystems A and B, one can trace out 
the subsystem B to obtain the reduced density matrix on 
subsystem A: pA ^ Trs \^) {^\- The Renyi entropies are 
defined as: 



Sn = z log (Trpl) 

1—71 

where n is an index parameter. Taking the limit n — > 
1; Sn recovers the definition of the usual von Neumann 
entropy. In this paper we will often discuss the Renyi 
entropy with index n ~ 2: 82 — — log {Tr (p^)) since it 
can be calculated most easily with the VMC methodic 
and at the same time, captures all the information that 
we are interested in. 

For a gapped phase in 2D with topological order and 
a disc shaped region A with smooth boundary of length 
La, the Area Law of the Renyi entropy gives: 

Sn = anLA - 7 (1) 

where we have omitted the sub-leading terms. Al- 
though the coefficient a„ of the leading 'boundary law' 
term is non-universal, the sub- leading constant 7, which 
is often dubbed as the TEE, is universal and a robust 
property of the phase of matter for which |$) is the 
ground state. When region A has a disc geometry, it 
has been shown that 7 for different degenerate ground 
states are identical and it is also insensitive to the Renyi 
entropy index ni^^. It equals 7 = log!?, where D is the 
total quantum dimension of the mode l^^'^'^ , and offers 
a partial characterization of the underlying topological 
order. 

However, when the subsystem A takes a non-trivial 
topology, or more precisely when the boundary of A is 
non-contractible, TEE contains more informationi^, as 
we will elaborate further in this paper. For simplicity of 
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illustration, throughout we focus on the case when the 
two-dimensional space is a torus T2 and the subsystem 
A wraps around the y direction of the torus and takes 
the geometry of a cylinder. For such a geometry, the 
?i'tli Rcnyi entropy corresponding to the wave function 
l*> Ej CjISj) is given by: Sn = UnLA-in-, where |Sj) 
is a special basis that we will describe in detail below and 

In 



In (fel) = 27 + ^ log \^p-df-^]^ (2) 

Here dj > 1 is the quantum dimension of the jth quasi- 
particle and pj = |cjp. For Abelian anyons, dj = 1. 
Note, dj shares the same subscript j as the states |Sj) 
because the states can be obtained by inserting a 
quasi-particle with quantum dimension dj (the ground 
state degeneracy on the torus is equal to the number of 
distinct quasi-particles) . This equation shows that the 
TEE for this geometry depends on the wave function 
through {pj} as well as the Renyi index n, unlike the 
case with disc geometry. 

What is the physical significance of the basis states 
I S j ) ? We claim that these are precisely the eigenstates of 
the nonlocal operators defined on the entanglement cut, 
which distinguish the topologically degenerate ground 
states. For example, in the case of the quantum HalU^ 
(Sec. IIIB 2p . these states are the eigenstates of the 
Wilson loop operator associated with the Chern-Simons 
gauge field around the hole exposed by the entangle- 
ment cut. Similarly, for a Z2 gauge theory (Sec. Ill CI 
), these are the states with definite electric and magnetic 
field fluxes perpendicular to the entanglement cut. For 
Abelian states, which have dj — 1 for all j and are the 
focus of this paper, the entanglement entropy associated 
with the states jSj) is minimum, i.e. heuristically, the 
entanglement cut has the maximum 'knowledge' about 
these states. For this reason we name them Minimum 
Entropy States (MESs). 



A. Strong Subadditivity and Topological 
Entanglement Entropy on the Torus: An 
'uncertainty' principle 

In this section wc discuss the TEE for bipartitions of 
a torus into two cylinders. This can be done by slic- 
ing the torus in two distinct ways, along the vertical or 
horizontal directions. Intuitively, one might expect both 
bipartitions would have the same TEE of 27, given the 
two disconnected boundaries of the cylinders. However, 
very general considerations based on strong subadditivity 
of von-Neumann entropy alone suggest that this expec- 
tation cannot be correct. In practice, it is known that 
for a wide class of topological phases, TEE of such non- 
trivial bipartitions indeed depends on the ground state 
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FIG. 1: A torus (the top and bottom sides and left and 
right sides are identified). Subregions A, B, C are defined 
as shown. Regions A and C are assumed to be well separated 
as compared to the correlation length. The regions AB and 
BC correspond to bipartitions of the torus into cylinders in 
orthogonal directions. 

selectedi^. Here we do not address ground-state depen- 
dence, rather we demonstrate that TEE cannot be iden- 
tical to its value for trivial bipartitions. It invokes strong 
subadditivity, a deep property of quantum information^^. 
This will allow us to come up with an uncertainty princi- 
ple, which constrains the amount of information we have 
when we cut the torus in two orthogonal directions. Its 
advantage is that it assumes almost nothing about the 
phase, except that it is gapped. 

Consider the ground-state wave function of a gapped 
phase in two dimensions and three non-overlapping sub- 
regions A, B, C. The von-Neumann entropies S follow 
the strong subadditivity condition^S: 

Sabc + Sb — Sab — Sec < (3) 

Note, this is only known to hold for von-Neumann 
entropies, not Renyi entropies in general. Now, con- 
sider a torus with subregions A, B, C as shown in the 
Fig. [T] Let us decompose the entropy into a part 
that arises from local contributions and a non-local TEE 
S = 5''°'=^' + S'*°P°. For a subregion with the topology 
of a disc, the TEE is expected to be 5*°p° = -7. Quite 
generally one can argue that 7 > utilizing the strong 
subadditivity conditioni^. For subregions defined on a 
simply connected surface, such as a disc, the TEE is pro- 
portional to the number of connected components of the 
boundary. If this was also true for the torus, we would 
expect 5*^^° = S^g^° = — 27 (since they have a a pair of 
boundaries). We now show this cannot be a consistent 
assignment of TEE on the torus. 

In order to isolate the topological part of the entropy, 
we assume that the regions A and C are well separated 
compared to the correlation length of the gapped ground 
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state. Then, the loeal eontributions cancel in the combi- 
nation above: S^""^ + 5']^'=^' - S^Xb"^ - S'J^g'' ^ 0. This 
can be argued following Refs.— i^^. For example, consider 
a local deformation near region A's boundary far away 
from the other regions. This change will be in S]^^q, but 
a nearly identical contribution will also appear in S^Xi3^7 
since it only differs by the addition of a distant region. 
These will cancel in the combination above. Thus, we 
can rewrite the Eqn. [3] as: 
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(4) 



This inequality implies the TEEs expected from the 
disc is not legal for the torus. 

For regions where the boundary is topologically trivial 
and contractible (such as ABC or B), one expects the 
TEE to be independent of the surface on which they are 



defined, and hence S" 
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AB and BC, whose boundaries wrap around the torus, 
are sensitive to the topology of the space they are defined 
on. Their TEEs satisfy: 



7BC + lAB < 27 



(5) 



where we have defined 8*^^°^^^ = —^ab{bc)- Clearly 
this does not allow both the TEEs to to be 27. In fact 
if one of them attains its maximal disc value, the other 
must vanish. Note, the TEE reduces the total entropy. 
Thus, when the entropy of a cut along one of the cycles 
of the torus attains its minimum value, i.e. we have most 
knowledge about the state on the cut, then along the 
orthogonal direction, the entropy associated with a cut 
must attain its maximal value, implying our knowledge 
is the least. Therefore this can be thought of as an uncer- 
tainty relation, between cuts that wrap around different 
directions of the torus. 



B. Ground State Dependence of TEE in a Chiral 
Spin Liquid 

In this subsection, to illustrate the state dependence 
of TEE, we study the entanglement properties in a lat- 
tice model of an SU (2) spin-symmetric CSL on a torus. 
The CSL has the same topological order as the half filled 
Landau level v = 1/2 Laughlin statei^ii^ of bosons (these 
bosons can be thought of as residing at the location of 
spin up moments), and has two-fold degenerate ground 
states on the torus. The topological order in CSL can 
be confirmed by calculating its topological entanglement 
entropy (TEE) numerically using Monte Carlo and veri- 
fying that it is non-zero and agrees with the field theoret- 
ical predictions^^. We note that here we are working with 
a generic wave function in this phase defined on a lattice, 
rather than with idealized zero correlation length states 
or the topological field theory. This will introduce new 



FIG. 2: The separation of the system into subsystem A, B, 
C and environment, periodic or antiperiodic boundary condi- 
tion is employed in both x and y directions, a: The subsys- 
tem ABC is an isolated square and the measured TEE has 
no ground state dependence, b: The subsystem ABC takes 
a non-trivial cylindrical geometry and wraps around the y 
direction, and TEE may possess ground stated dependence. 



conceptual issues - in particular, the connection between 
MESs and lattice ground states will be discussed. 

We begin by reporting the results of a numerical ex- 
periment. We extract TEE of linear combinations of the 
two ground states of the CSL, and show that it indeed 
depends systematically on the chosen linear combination, 
when the entanglement cut wraps around the torus. We 
will then predict theoretically the dependence and find 
excellent agreement as shown in Fig. [31 



1. Numerical Study of Ground State Dependence of TEE in 
Gutzwiller Projected CSL states 

Wave functions of an SU(2) spin symmetric CSL are 
obtained in the slave particle construction. We write 
the spins as bilinear in fermions S = \fl[(T]aa' fa' 
and assume a chiral d-wave state for the fermions. 
Operationally, the spin wavefunctions are obtained by 
Gutzwiller projection of a d^^-yi + idxy superconductor 
to one fermion per site. More technical details regard- 
ing this wave function are in Appendix [BJ We consider 
the system on a torus. Before projection, one can write 
down different fermion states, by choosing periodic or 
anti-periodic boundary conditions along x and y direc- 
tions. These boundary conditions are invisible to the 
spin degrees of freedom which are bilinear in the fermions 
and lead to degenerate ground states^. We denote the 
ground states by the mean field fluxes in x and y direc- 
tions as \lpi, tp2), (pi,2 = 0, TT. The two fold degeneracy of 
the CSL implies that only two of the four ground states 
|0,0), |7r,0), |0,7r), |7r,7r) are linearly independent. Here 
we consider linear combinations of |0, tt) and \-k, 0), which 
we have numerically checked to be indeed orthogonal for 
the system sizes that we consider: 



= cos(/i|0,7r) -hsin0|7r,O) 



(6) 
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We calculated TEE for the state |$) using VMC 
method and Gutzwiller projected wave functions based 
on Eqn IBll An efficient VMC algorithm which allows to 
study a linear combination of Gutzwiller projected wave 
functions was developed and detailed in Appendix [Xj To 
our knowledge, this is the first numerical study to accom- 
plish this. 

The geometry and partition of the system arc shown 
in Fig. [2]3. The total system size is 12 lattice spacings 
in both directions with rectangles A and B being 6x4 
and rectangle C f 2 x 4. Note that the subsystems AC, 
BC, AB, C and ABC all wrap around y direction so that 
their TEE will all be equal (and denoted 7'). This is the 
quantity we wish to access. For contractible subsystems 
A and B it remains the same as that expected for a region 
with a single boundary, cut out of a topologically trivial 
surface (such as a bigger disc) 7. Wc use the construction 
due to Kitaev and PreskilU^ and effectively isolate the 
topological contributions in the limit of small correlation 
length, by evaluating the combination of entropies Sa + 
Sb + Sc — Sab — Sac — Sbc + Sabc- This combination 
is related to the TEE by: 

-27 + 7' = Sa + Sb + Sc 

— Sab — Sac — Sbc + Sabc 

= 25"^ — 2Sac + Sabc (7) 

In the second line we have exploited symmetries of 
the construction to reduce the problem to calculation of 
the Rcnyi entropy 6*2 of three regions A, AC and ABC 
for each (f). To measure S2 numerically, we calculated 
the expectation value of a SwapA operator, see Refi^ 
for an elaboration of the method used. Our results for 
27 — 7' ((/)) corresponding to different linear combinations 
parameterized by 4> are shown in Figl3] This is one of the 
main results of this work. 

We note the TEE strongly depends on the particular 
linear combination chosen. The zero of the curve implies 
that the TEE 7' = 27, intuitive value for an entangle- 
ment cut with two boundaries. The corresponding state 
is the MES. We note that the MES occurs at a nontrivial 
angle. Understanding this requires connecting the lat- 
tice states and the field theory which is done below. Wc 
predict this angle to be 0.1257r and the overall TEE de- 
pendence to be Eqn. [Til which is plotted as the solid 
curve Fig. [3l in rather good agreement with the numeri- 
cal data. 



2. Theoretical Evaluation of Ground State Dependence of 
TEE m CSL wavefunctions 

A calculation of ground-state dependence of TEE in- 
volves two steps. In the first step, we ask the following 
question: given a state expressed as a linear combination 
of MESs, what is the expected TEE? For the CSL, this 
question has already been answered by Ref. [ij that TEE 
for a state = ai|l) -I- a2|2) is 



0.8- 




FIG. 3: Numerically measured TEE 27 — 7' for a CSL ground 
state from linear combination 1$) = cos0|O, tt) + sm<j}\n,0) 
as a function of <j) with VMC simulations using geometry in 
Fig. ^jp. The solid curve is the theoretical value from Eqn. 
111! The periodicity is 7r/2. 



7'-27 = log(|ai|4 + |a2|4) (8) 

where |1), |2) are MESs for cutting the torus in the 
direction in question. 

Second, wc need to understand the relation between 
the MES and the physical states that appear in the 
Gutzwiller wave function. In general it appears that the 
only way to identify MESs in a generic wave function is 
by calculating the TEE. However, when the lattice model 
has additional symmetry, that can also be used to iden- 
tify MESs. Here, we have a 12 x 12 system defined on a 
square lattice and we will exploit the tt/2 rotation sym- 
metry to establish a connection between the flux states 
\ipi, 1P2) of the Gutzwiller ansatz and the MESs. 

The Gutzwiller projected ground states of the CSL, 
|0, 0) and |7r,7r) are clearly invariant under a 7r/2 rota- 
tion symmetry upto a phase factor. A simple calculation 
shows that the |0, 0) state acquires phase factor —1 while 
the Itt, tt) state acquires no phase under rotation. Sim- 
ilarly, the states (|0,7r) ± |7r,0)) acquire a phase ±1 
under rotation. Having established the transformation 
of lattice states under rotation, we now study how the 
MESs in the field theory respond to rotations. Wc will 
see that tt/2 rotation in the basis of the MESs is de- 
scribed by the modular S matrix. The eigenvectors of 
the modular S matrix will then be identified with lattice 
states that are rotation eigenstates. 

The CSL has the same topological order as the half 
filled Landau level = 1/2 Laughlin stat e^^'^^ of bosons. 
The field theory describing the topological order oi a. v = 
1/k Laughlin state is described by the following Chern- 
Simons action. Note, here only the very long wavelength 
degrees of freedom are retained: 
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One can define the Wilson loop operators Ti = e = 
^ifa^dx g^j^j rj.^ ^ giSa ^ g^fa^rfy around the two distinct 

cycles of the torus. In terms of 9i, the action is given by 



S = i— I (116162 
27r 



which implies that at the operator level [01,^2] = 



or 




T1T2 = TaTie^''*/^ 

Owing to the above relation, there are k orthogonal 
ground states IV'm) that can be chosen to transform under 
Ti as 



„2iri(m-l)/fc 



l^m+l) 

In the case of a CSL phase, k = 2. Let us label the two 
3generate groum 
eigenstates of T2 : 



degenerate ground states as (1, 0)^ and (0, 1)^, which are 



FIG. 4: The presence of two fluxes 9i and 62 on a torus. In 
cylindrical cut A that wraps around y direction, only 62 is a 
measurable. 



deriving S matrix for CSL state using the entanglement 
properties of the ground states in Sec. IIII Al Here we 
restrict ourselves to the calculation of TEE for the CSL. 

Since we are interested in the entanglement entropy 
with respect to the cut shown in Fig. 21 let us represent 
all our states in the basis of the eigenstates of T2, i.e. the 
states (0, 1) and (1, 0). Then, by matching eigenstates of 
the S matrix in the above basis and rotation eigenstates 
of the lattice problem, we conclude: 



72(1, of = (l,Of 

r2(o,if = -(o,if 
ri(i,of = (o,if 
ri(o,if = (i,of 

The last two equations are due to the commutation 
relation T1T2 = —T2T1. It follows that the eigenstates of 
Ti are (1, if /V2 and (1, -1)^ /V2. 

The significance of the Ti^2 eigenstates is that they are 
MESsi^, for cuts whose boundaries are parallel to the 
loops used to define Ti^2- This is because eigenstates of 
these loop operators have a fixed value of flux enclosed 
within the relevant cycle of the torus, which minimizes 
the entanglement entropy for a parallel cut as shown in 

Fig. m 

Now consider a 7r/2 rotation, under which 61 — >■ 62 and 
62 -61 so Ti T2 and T2 Tf ^ = Ti. Thus, the 
matrix representing the effect of 7r/2 rotation for CSL in 
T2 cigenstatc basis is: 

S=(^ (9) 

V V2 J 

Note that we have used the symbol S for the above 
matrix because it is indeed the modular S matrix of 
the Chern-Simons topological quantum field theory cor- 
responding to a CSL. We recall that the modular S ma- 
trix transforms the eigenstates of one Wilson loop opera- 
tor T2 to those of Ti . We will return to the discussion of 




We can now expand the general linear combination 
state in MESs: 

1$) = cos(/)|0,7r) +sin0|7r,O) 

= (cos(0-|),sin(0-^)) (10) 

Then, according to Eqn. ^ , theoretically one expects 
the following expression for TEE: 



27-7' ^ log . (11) 
3 + sm (40) 

which is compared with the numerical data in Fig. |3l 
The MES occur at the value of 7r/8 (mod ■n/2). 



C. Toric code and Z2 spin liquid 

In this subsection we use the Kitaev's Toric code 
model^ as a pedagogical example to understand ground- 
state dependence of TEE and the nature of the MESs for 
a Z2 gauge theory. 

Consider the toric code Hamiltonian of spins defined 
on the links of a square lattice^: 



p 



FIG. 5: Illustration of a lattice of the toric code model, the 
links spanned by star and plaquette are highlighted in red and 
blue, respectively. 



H 



(12) 




FIG. 6: A snapshot of the ground state on the cylinder. 
Closed-loop strings ( "Z2 electric fields" ) can wrap around the 
cylinder. The ground states are doubly degenerate, corre- 
sponding to even and odd winding number sectors. The to- 
tal number of string crossings the cut A equals the winding 
number, modulo 2. The number of string crossings at the 
boundary F is even in the degenerate ground states. 



where s and p represent the links spanned by star and 



plaquette as shown in the FiglH and As 



Bp = rijGp "^j ■ Since all individual terms in the Haniil- 
tonian commute with each other, ground states are con- 
structed from the simultaneous eigenstates of all As and 
Bp. Define the operator W^{C) associated with a set of 
closed curves C on the bonds of the lattice, as follows 



VF^(C) 



(13) 



Then the ground state is an equal superposition of all 
possible loop configurations: J2c^ab(.^)\^^^x), where 
Ivac^,) is a state with = —1 on every site. The closed 
loops are interpreted as electric field lines of the Z2 gauge 
theory. We now consider two geometries, first the cylin- 
der and then the torus. The former case has a pair of 
degenerate ground states, and is the simplest setting to 
demonstrate state dependence of TEE. 



Cylinder Geometry 



/^E(K},o)i*f.},o) 
7^E(K},o>i*f.},i> 



I* 



B 



{*},l)l*{<?,},l 



{qi} 



(14) 



where the Nq distinct configurations represented by 
{qi} denotes the electric field configurations at the cut. 
The number of field lines crossing the cut is always even, 
since the ground state is composed of closed loops. For 
trivial bipartitions, this exhausts all terms in the Schmidt 
decomposition^^. However, given that the boundary of 
the cut is non-contractiblc, the additional index 0, 1 ap- 
pears which counts the parity of electric field winding 
around the cylinder, within a partition. These are cor- 
related between the two partitions, for the fixed winding 
number ground states. This is the key difference from 
a trivial bipartition, leading to the ground-state depen- 
dence of TEE. 

We now calculate the entanglement entropy associated 
with such a cut for an arbitrary linear combination of 
these two ground states IvP) = co|^o) + ci|^i), with unit 
norm. Using Eqn. [14] one can easily verify: 



On a cylinder, the Hamiltonian in EqnlT^] leads to 
a pair of degenerate ground states (the As part of the 
Hamiltonian is suitably modified at the boundary of the 
cylinder to only include three links). The two normal- 
ized ground states |^o), ^-re given by equal superpo- 
sitions of electric field loop configurations which have an 
even and odd winding number around the cylinder re- 
spectively (see Fig. [6]). Consider now partitioning the 
cylinder into two cylindrical regions A and B. Then 
the Schmidt decomposition of these ground states can 
be written as: 




+ (co-cOlvpf^^j^Jlvff^^j,.)] (15) 

-here |*f f^^^) ^ (l*^*^)^) ± /V2. For a 

Schmidt decomposition l^-) = J2a \/^l*a)l*f ) the nth 
Renyi entropy is given by: 5„ = j^Jog (Y^a^a)- We 
arrive at: Sn = log A^^~"[p" + pi], where p± = 
\co ± cip/2. Recognizing that the closed loop constraint 
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leads to Nq = 2^ ^, where L is the length of the cut, and 
using the definition of TEE in Eqn. [T] we have: 



7;=log2-^log(p!^+p!^) (16) 

Thus, for the electric field winding eigenstates |Co.i) 
where p± = 1/2, the TEE vanishes. However, for their 
equal superpositions when one of p+ or p_ vanishes, the 
TEE attains its maximal value log 2. These are eigen- 
states of the Wilson loop operator that encircles the 
cylinder and measures the Z2 magnetic flux (vison num- 
ber) threading it. An example of such a flux operator is 
F = rijgQ '^h where Q is a closed curve that loops once 
around the cylinder, such as the boundary T in Fig. 1^1 
Since TEE reduces the entanglement entropy, the maxi- 
mum TEE states correspond to MESs. Why these MESs 
are eigenstates of flux through the cylinder for this par- 
ticular cut? The number of electric field lines crossing the 
boundary F is always even. This constraint carries some 
information and hence lowers the entropy by bringing in 
the standard TEE of log 2. On the other hand, the topol- 
ogy of the cut boundary F allows for a determination of 
which magnetic flux sector the cylinder is in. A state 
that is not an eigenstate of magnetic flux through the 
cylinder leads to a loss of information and hence a posi- 
tive contribution to the total entanglement entropy (and 
reduces TEE). This suggests that the MESs are eigen- 
states of loop operators which can be defined parallel to 
the cut F. This is further substantiated by the result for 
the torus case discussed below, where they are simulta- 
neous eigenstates of magnetic flux enclosed by the cut 
and electric flux penetrating the cut. 



2. Torus Geometry 

The four degenerate ground states are distinguished 
by the even-odd parity of the winding number of electric 
field lines around the two cycles of the torus. The oper- 
ator W^{C), which generates the set of closed loops C 
can be used to write the ground states: 

c 

where the subscript a (b) takes on binary values 0, 1 
and denotes whether the loops C belong to the even or 
odd winding number sectors along the x (y) direction, 
and |vac^) is a state with = —1 on every site. The 
four ground states cannot be mixed by any local operator 
and hence realize a Z2 topological order. Let us consider 
a ground state as the following linear combination: 

a,b=0,l 



We are interested in calculating entanglement entropy 
for the state |^) corresponding to the partition shown in 
the FigHland the dependence of TEE on parameters Ca,b- 
After straightforward algebra (see details in Appendix 
[C|) . one finds the following expression for subsystem A 
with boundaries of length L: 

Sn = ilog(2) - 7; 

where 



7; = 21og(2)-^log5]p," (18) 

and 



Coo 


+ coi 






2 




|coo 


- Coi 






2 




|cio 


+ Cll 






2 




|cio 


- Cll 





(19) 



This is indeed consistent with Eqnl^l given that 7 = 
log D = log 2 and dj = 1 for an Abelian topological order 
with = 4 degenerate ground states. Further, EqnlTSl 
readily leads to the following four MESs: 

iSi) = -i=(ieoo) + ieoi)) 
= -i=(ieoo) - icoi)) 

|S3)-^(|ClO> + |Cll)) 

is4) = -i=(ieio>-iai» (20) 

What is the physical significance of these four states 
being the MESs? Similar to the cylinder geometry case, 
these states are the simultaneous eigenstates of Wilson 
loop operators that encircles the torus and measures the 
Z2 magnetic and electric fluxes threading it, as shown 
in Table H] and Fig. We leave more detailed algebra to 
Appendix [Cl 

When 7^ is maximized, the corresponding Sn is min- 
imized, providing the maximum possible information 
about a given state. Since the cut is made along y, it 
can measure the Z2 magnetic and electric fluxes directed 
parallel to x. Hence the MES \Ea) with definite mag- 
netic and electric flux sectors, maximizes the TEE with 
'Jtopo = 2 log(2), a contribution of log(2) from each of the 
two boundaries. As one linear superposes different MESs 
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MES 


Ty 


Fy 


quasi-particle 










1 


H2 


1 





m 


23 





1 


e 


S4 


1 


1 


em 



TABLE I: List of Z2 magnetic flux Ty, Z2 electric flux Fy and 
corresponding quasi-particle of Wilson loop operator for the 
four MESs I Ha) of the toric code with system geometry in 
FigSl The definitions of Ty and Fy are in Appendix [Cl 

|So,), the mformation obtained from measuring magnetic 
and electric sectors becomes scrambled; especially, at the 
extreme case of equal superposition of |Sa), all informa- 
tion about the global quantum numbers has been lost 
and wc have 7' = 0. This offers another example where 
MESs are the eigenstatcs of loop operators defined on the 
cylinder from the entanglement cut. 

Numerical Verification in a Generic Z2 wave- 
function: To provide an example of the state depen- 
dence of TEE in a Z2 topologically ordered beyond the 
simple toric code model, we employ our VMC method to 
calculate the TEE of an SU{2) symmetric Z2 spin liq- 
uid for a cylindrical bipartition of the torus. This is a 
generic spin wave function, and has a finite correlation 
length unlike the toric code. It is obtained by Gutzwiller 
projecting a mean-field BCS state on a square lattice (see 
Appendix [B] for details). 

For the Z2 spin liquid wave function obtained by 
Gutzwiller projecting the above mean-field ground state, 
a calculation of TEE corresponding to a disc shaped re- 
gion A was performed in Refji^. It was established that 
the state is indeed topologically ordered with the value 
of 7 sa 0.584 ± 0.089, close to the expected theoretical 
value of 7 = log 2 k, 0.693. Here we are concerned with 
the calculation of TEE for a cylindrical bipartition. We 
study systems with 12 lattice spacings in the x direction 
and 8 lattice spacings in the y direction. The subsystem 
separation scheme is similar to Figl5}D, where A and B 
are 4x4 squares and C is a 4 x 8 rectangle, see Fig. [T] 
The TEE 27 — 7' can be evaluated in a similar way as 
EqnEl 

By employing periodic boundary conditions in both 
directions, the state we study has zero magnetic fiux 
along both directions, but it is not a MES. Instead, it 
is an equal superposition of two MESs states that have 
zero magnetic flux and or 1 electric flux respectively 
through the cylinder. Therefore, according to EqnfTSl 
with pi = P3 = 1/2, p2 = Pi = and Rcnyi index n = 2, 
one finds: 

4 

7' - 21og(2)+log(^p2)^lQg(2) 

27-7' = log (2) 

The VMC simulation yields 27 - 7' = 0.576 ± 0.108, 
close to the expected theoretical value log(2) ^ 0.693. 
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1 













FIG. 7: The separation of the system into subsystem A, B, 
C and environment, periodic boundary condition is employed 
in both X and y directions. The subsystem ABC takes a non- 
trivial cylindrical geometry and wraps around the y direction. 

The smaller than expected values for both '^ 1^^'^^ and 7' 
are probably due to quasi-particle excitations across a fi- 
nite gap, causing breaking of Z2 electric field lines over 
the finite system size we consider. Indeed, spin correla- 
tions decay slower for the Z2 state, as compared to the 
CSL. which also arrives closer to its expected value of 
TEE for both cylindrical bipartition as well as for the 
disc shaped region A^. 



III. EXTRACTING STATISTICS FROM 
TOPOLOGICAL ENTANGLEMENT ENTROPY 

The modular S and lA matrices describe the action 
of certain modular transformations on the degenerate 
ground states of the topological quantum field theory. 
On the other hand, the braiding and statistics of quasi- 
particles are encoded in the S and lA matrices. For 
Abclian phases, the ij'th entry of the S matrix corre- 
sponds to the phase the «'th quasi-particle acquires when 
it encircles the j'th quasi-particle. The U matrix is di- 
agonal and the w'th entry corresponds to the phase the 
I'th quasi-particle acquires when it is exchanged with an 
identical one. Since the MESs are the eigenstates of the 
nonlocal operators defined on the entanglement cut, the 
MESs are the canonical basis for defining S and U. The 
modular matrices are just certain unitary transforma- 
tions of the MES basis. As argued in Appendix |D1 the 
S matrix acts on MESs as an operator that implements 
7r/2 rotation while the US matrix corresponds to 27r/3 
rotation of MESs. 



A. Modular S-matrix of CSL from TEE 

Let's consider the CSL wave functions studied in Sec. 
Ill B 11 and assume that we did not have any information 
about the individual quantum dimensions or the mod- 
ular S matrix. The only information that is provided 
is the two-fold degenerate ground-state wave functions 
Itt, 0) and |0,7r). We construct the linear combination 
1$) as EqnfTUland calculate its TEE for a non-trivial bi- 
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partition as Fig|4] on a 7r/2 rotation symmetric lattice. 
Consequently, we get the 27 — 7' dependence on param- 
eter (j) in FiglSl 

Wc notice that the minimum of the 27' — 7 attained 
is approximately zero. According to Eqnl^l this implies 
that at least one of the quantum dimensions di should be 
1. Since the total quantum dimension D 

^72, this implies that 
the MES lies at 



dl 



- d"^ 
"1/2 



di/2 = 1- Also, we see that 
0.147r by fitting FigOlto EqnH 
For system with square geometry the S matrix de- 
scribes the action of 7r/2 rotation on the MESs. Since 
the two states |0,7r) and |7r,0) transform into each 
other under 7r/2 rotation, this implies that in the basis 
{|0, tt), Itt, 0)}, the modular S matrix is given by the Pauli 
matrix cr^. To change the basis to MESs, we just need 
a unitary transformation V that rotates |0, tt) Itt, 0) basis 
to the MESs basis. The V is determined by the fact that 
one needs to rotate |0, tt) [tt, 0) basis by an angle « 0.147r 
to obtain MES (this is the numerically determined value, 
the exact value being tt/8). Therefore, 



error in pinpointing the MES and thereby leading to a 
more accurate value of the S matrix. As another applica- 
tion, we study the action of modular transformation on 
the MESs ISq) for the Z2 gauge theory in Appendix [El 



B. Algorithm for extracting modular S matrix 
from TEE 

In the last subsection and Appendix [E] we calculated 
the modular S matrices for the CSL and Z2 toric code 
model respectively using the transformation properties of 
the basis states \^ab) under 7r/2 rotation transformations 
TZ and translating it into the canonical basis of MESs 
I Si): S = WTZU. However, such tt/2 rotation symmetry 
is not necessary. Even without symmetry, it is possible 
to obtain the S matrix by studying the modular transfor- 
mation between certain nontrivial pair of sets of MESs. 

Starting from the definition, the modular S matrix has 
the following expression: 



5 = 



1 

1 



V 



where 



V: 



cos(0.147r) 
sin(0.147r) 



-sin(0.147r)e*'^ 
cos(0.147r)e''^ 



from the two MESs: (cos(0.147r), sin(0.147r))'^ and 
(- sin(0.147r),cos(0.147r))'^ and tp is an undecided phase. 
This yields the following value for the approximate mod- 
ular S matrix 



sin(0.287r) cos(0.287r)e*'^ 
cos(0.287r)e-*'^ - sin(0.287r) 



The existence of an identity particle requires positive 
real entries in the first row and column and implies ip — 0, 
which gives: 



0.77 0.63 
0.63 -0.77 



Comparing this result with the exact expression in 
Eqn. ini wc observe that even though the S matrix ob- 
tained using our method is approximate, some of the 
more important statistics can be extracted rather exactly. 
The above S matrix tells us that the quasi-particle cor- 
responding to do = 1 does not acquire any phase when 
it goes around any other particle and corresponds to an 
identity particle as expected, while the quasi-particle cor- 
responding to di/2 = 1 has semion statistics since it ac- 
quires a phase of tt when it encircles another identical 
particle. Numerical improvements can further reduce the 



5, 



Q,3 



D 



7^ \7:V 



(21) 



Here D is the total quantum dimension and x and y 
are two directions on a torus. Ean l21l is just a unitary 
transformation between the particle states along differ- 
ent directions. In the case of a square system as in the 
last subsection, the S matrix acts as a 7r/2 rotation on 

the MES basis In general, however, x and y do 

not need to be geometrically orthogonal, and the system 
docs not need to be rotationally symmetric, as long as 

the loops defining |Sj|j) and interwind with each 

other. Therefore, the modular S matrix can be derived 
even without any presumed symmetry of the given wave 
functions. Note that there is an undetermined phase for 

each and 

the rows (columns) , which may be fixed by the existence 
of an identity particle. 

Let's start with the two primitive vectors wi and 1^2 
and determine the transformation of the MESs of W2 to 
those of w'2 given by: 



) , therefore a phase freedom between 



niWi 
n2Wi 



miW2 

m2W2 



(22) 



With nim2 — min2 = 1 by definition of the modular 
transformation. We restrict 712 — —1, which means the 
cross product: 



W2 X W2 



-W2 X U'l ~ Wi X W2 



A 



(23) 



A is the surface area of the torus. 

The corresponding modular matrix can be expanded 
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f ni 1~ 1211712 \ ^ f ^ \ ( ^ \ ( ^ ^"^2 ^ 

\-l m, ) \ 1 )\ I ) 

^ u-"'su-"'' 

Correspondingly, according to Appendix |D] the trans- 
formation: 

Because U matrix is diagonal by definition, its left 
(right) matrix products only adds an additional phase 
factor to each row (column) and can be eliminated. 
Therefore, without any argument on the symmetry, the 
generalized algorithm: 

1. Given a set of ground state wave functions |^q), 
calculate the TEE of an entanglement bipartition 
along W2 direction, for a linear combination |$) = 
'^Cal^a) ■ Search for the minimum of TEE 27 — 7' 
in the Cq parameter space. That gives one MES 
1 5^) and the corresponding quantum dimension 
21og((i^) ~ 27 — 7'. Note that the existence of 
an identity particle ensures at least one minimum 
TEE 27 - 7' = 0. 

2. iterate step 1 but with Ca in the parameter space or- 
thogonal to all previous obtained MESs Con- 
tinue this process until we have the expressions 
for all IS/3). This gives a unitary transformation 
matrix Ui with the a/3'th entry being Cq^, which 
changes the basis from \^a) to IS/3). Note that there 
is a relative U{1) phase degree of freedom for each 
IS,). 

3. Repeat step 1 and step 2 but with the entanglement 
cut along direction, which satisfies Eqn. 1221 and 
Eqn. 1231 and obtains the unitary transformation 
matrix U2- 

4. The modular S matrix is given by U2^Ui except for 
an undetermined phase for each MES correspond- 
ing to a row or a column. The existence of an iden- 
tity particle that obtains trivial phase encircling 
any quasi-particle helps to fix the relative phase be- 
tween different MESs, requiring the entries of the 
first row and column to be real and positive. This 
completely defines the modular S matrix. 

The above algorithm is able to extract the modular 
transformation matrix S and hence braiding and mu- 
tual statistics of quasi-particle excitations just using the 
ground-state wave functions as an input. Further, there 
is no loss of generality for non-Abelian phases, which can 
be dealt by enforcing the orthogonality condition in step 
2 which guarantees that one obtains states with quantum 
dimensions da in an increasing order. 

In Appendix |E] we take the square lattice toric code 
model as an example once again, but without presuming 
any symmetry of the system. 



C. Extracting other modular matrices from TEE 

In Appendix |E1 we calculate the U matrix for the Z2 
toric code model, given the simple action of W on l^ob)- 
Though we were unable to find a general algorithm for 
the U matrix, as we did for the the S matrix in the last 
subsection, in the presence of certain symmetries U can 
indeed be extracted given a set of ground-state wave func- 
tions I^q). This is achieved by first calculating the action 
TZ on the states |^q,) under this symmetry operation, and 
then translating it into the action on MESs. Specifically, 
the corresponding modular matrix is given by U^TZU, 
where the unitary matrix U is obtained through the first 
two steps of the algorithm in the last subsection. 

The aforementioned symmetry to extract S matrix is 
the 7r/2 rotation, as shown in Sec. IIII Al and the first 
example in Appendix [Ej but it may be generalized to 
symmetries such as rotation of other angles and even re- 
flection symmetry (see Appendix |D]) . More interestingly, 
when the symmetry operation 7?, is a 27r/3 rotation, one 
gets the US matrix. Hence, if one starts with an arbi- 
trary basis for the degenerate ground state manifold 
of a topological order, the problem of S and U matrices 
can be reduced to the transformation property of cho- 
sen basis states |^q) under tt/2 and 27r/3 rotations and 
the unitary transformation that translates |^q,) basis to 
the MESs |Sq). To illustrate this point, we extract the 
US matrix for the Z2 gauge theory in Appendix [E] by 
putting the Z2 toric code on triangular lattice which has 
27r/3 rotation symmetry. 

IV. CONCLUSION 

In this paper, we studied two topologically ordered 
phases, the chiral spin liquid (CSL) and the Z2 spin liq- 
uid using VMC method for Gutzwiller projected wave 
functions numerically and the Z2 toric code model ana- 
lytically, and showed that the topological entanglement 
entropy (TEE) depends on the chosen ground state when 
the entanglement bipartition is topologically nontrivial. 
We also determined the minimum entropy states (MESs) 
and explained their physical significance. 

As an application of the physical significance of MESs, 
we suggested an algorithm for extracting the modu- 
lar transformation matrices that determine the statistics 
such as quasi-particle self and mutual statistics and in- 
dividual quantum dimensions. These matrices also de- 
termine the central charge of the edge state modulo 8^^. 
Given the ground states, this algorithm determines the 
topological order to a large extent. We note that Wen 
proposed a different way to extract S and U matrices by 
calculating the non-abelian Berry phased which in prac- 
tice may be difficult to implement, especially on a lattice, 
since it requires calculating the degenerate ground states 
ipn of the system as function of the modular parame- 
ter r ~ LO2/UJ1 and calculating the derivatives such as 
(V„(r)|9,|^™(T)). 
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We note that there may be cases where the 7r/2 and 
27r/3 rotations of the MESs may not be exactly iden- 
tifiable with the modular S and L{S matrices respec- 
tively. This may happen, for example, when the particles 
have an internal angular momentum which may cause 
the wavefunction to acquire an additional phase upon 
rotation, over and above the phase due to underlying 
topological structure. If the MESs correspond to spin- 
singlet spin-liquid wave functions (such as CSL studied 
in this paper) and/or string- net models (such as toric 
code model) where there is no such internal structure, 
there should not be such additional phase. Further, since 
all MESs are locally same, they should all acquire same 
extra phase due to any local physics and therefore, the 
extra phase may be separable from the topological phase. 

For quantum Hall systems, because of the bulk-edge 
correspondence, the fusion algebra and topological spin 
of the bulk quasi-particlcs also determine the fusion rules 
and scaling dimensions for the primary fields in the chiral 
CFT at the edge. Therefore, in the context of quantum 
Hall systems, the entanglement entropy of the ground 
state manifold determines robust features of the fields in 
the corresponding edge CFT. 

It would also be interesting to consider generalization 
of the methods developed in our paper to higher dimen- 
sions. Discrete gauge theories furnish the best known 
theories with long-range entanglement in Z? > 3 dimen- 
sions and akin to D = 2, they again support degenerate 
ground states on the torus. In a recent paper—, it was 
shown that these theories again have non-zero TEE that 
is proportional to log(|G|), the number of elements in 
the gauge group. A simple generalization of the method 
developed in this paper shows that TEE for bipartition 
that has non-contractible boundaries will again depend 
on the ground state, and one will again find certain MESs 
that have the maximum knowledge of the quantum num- 
bers associated with an entanglement cut. Yet we are 
not aware of simple generalization of modular transfor- 
mations to higher dimensions, the meaning of the matrix 
that relates MESs for orthogonal entanglement cuts in 
higher dimensions requires further investigation. 
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Appendix A: Variational Monte Carlo method for a 
linear combination of wave-functions 

To calculate TEE for wave functions of different linear 
combinations, it is important to establish a VMC algo- 
rithm for wave function as 1$) = cos(j)\^i) +sin(/)|$2), 



where we assume |$i) and |$2) are properly normalized. 
In our case, |<I>i) and |$2) are two degenerate ground 
states, (a|$i) and (a|<I>2) are single Slater determinants 
products for each configuration \a), making (q;|$) a sum 
of two Slater determinants products. However, it may 
also be generalized to the situation of any wave func- 
tions. 

In the VMC scenario, the central quantity to evaluate 
in each Monte Carlo step is the ratio of (a'|<I>) / (a|<I'), 
which now has the form: 

(a'|$) _ cosi?;)(q;'|$i) +sin(?;)(a^|$2) 
(a|$) cos</i + sin(/> (a|$2) 

It is usually much less costly to calculate ratio of 
(a'l^i) / (a|$i) and (a'|$2) /(a|^2) if la) and \a') are 
locally different. For our case, when \a) and \a') differ 
only by one spin(electron) exchange, a much less costly 
and more accurate algorithm may be implemented for 
the ratio of determinants with only one different row or 
column. Unfortunately, after linear superposing different 
!<!>;), Eqn IAlI no longer has such a privilege. 

However, one can re express Egn lAll as: 

(a'|<i>) a + 6c-tan(/) 

(q|<I>) l + c-tan(/!) 

where 

a - (a'|$i)/(a|$i) 

b = (a'|$2)/(a|$2) 

are again ratio of determinants and can be effectively 
evaluated, and 

c= (a|$2) / (a|$i) 

can be efficiently kept track of with c' = a^^hc when- 
ever the update |a) — > \a') is accepted in a Monte Carlo 
step. In practice, numerical check should be included to 
make sure error for c does not accumulate too much after 
a certain number of Monte Carlo steps. 

This algorithm may be easily generalized to the linear 
combination of n wave functions, with the computational 
cost only n times that for a single wave functions. 

Appendix B: Variational wave functions for Chiral 
Spin Liquid and Z2 Spin Liquid. 

Chiral Spin Liquid Prom Gutzwiller Projection: 

The lattice wave function for the CSL states that we con- 
sider are obtained using the slave-particle formalism by 
Gutzwiller projecting a d~\-id BCS statei^i^. Specifically, 
we Gutzwiller project the ground state of the following 
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FIG. 8: Illustration of a square lattice hopping model con- 
nected with a d + id superconductor. While the nearest 
neighbor hopping is along the square edges with amplitude t 
{{—t) for hopping along dash lines), the second nearest neigh- 
bor hopping is along the square diagonal (arrows in bold), 
with amplitude +iA (— iA) when hopping direction is along 
(against) the arrow. The two sublattices in the unit cell are 
marked as A and B. 



This mean field model is readily solvable, with disper- 
sion: 



Ek = 



|A|| 



Afc 



2 (cos [kx] + cos iky)) 
2?] [cos {kx + ky) + cos (fc; 



-2iA [cos {kx + ky 



ky)] + ( 

COS {kx ~ ky)] 



We choose = A = 1.5 for a large gap and our cal- 
culation estimates that the correlation length is as short 
as ^ 1.3 lattice spacings. Please refer to Refi^ for more 
details on the construction. 



Hamiltonian of electrons hopping on a square lattice at 
half filling: 



H = Y^t 

{ij 



ikc]ck 



(Bl) 



Here i and j are nearest neighbors and the hopping 
amplitude tij is t along the ij direction and alternating 
between t and — i in the x direction from row to row; and 
i and k are second nearest neighbors connected by hop- 
pings along the square lattice diagonals, with amplitude 
iAifc = iA along the arrows and iA^fe = — iA against the 
arrows, sec FiglS) The unit cell contains two sublattices 
A and B. This model leads to a gapped state at half filling 
and the resulting valence band has unit Chern number. 
This hopping model is equivalent to a d -I- zd BCS state 
by an SU{2) Gauge transformation. We take A = 0.5t 
to maximize the relative size of the gap and minimize the 
finite size effect. Please refer to Re£i^ for further details 
regarding the exact form of the wave function. 

Z2 Spin Liquid from the Gutzwiller Construc- 
tion: 



H 



(ij) 



h.c.) + ^^|4tV» 



where ip, = (^f^, /| j 



^1,2,3 g^j.g pg^nii matrices. The 



second term is related to chemical potentials, we set 
a^'^ = 0, with aj fixed by the conditions (Jt/'V^'^''^?/') = 0. 
Matrices fiij connecting nearest and next nearest neigh- 
bors: 



Appendix C: Minimum entropy states of Toric code 
model on dividing torus 



In this appendix we Schmidt-decompose the individual 
Toric code ground states \'^) in Eqn[T7]for the bipartition 
of a torus in FigUl It is helpful to introduce a virtual cut 
A which wraps around the torus in the x direction, and 



defir 



5- 



MB) 



{11} b/ normalized equal superposition 

of all the possible configurations of closed-loop strings 
C in the subsystem A (B) with the partition boundary 
condition specified by {qi = 0,1}, I = l,2,...,i(so L is 
the total length of the boundary), and the number of 
crossings of the virtual cut A modulo 2 equals b ~ 0, 1. 
The four ground states may now be expanded as 



1* 



« {9i}ea 
I* 



{gj,0>l*{9,},(, 



{gi},l/l*{g,},(6+l)mod2/ 



Here {qi} & a = (1) denotes that the only even (odd) 
number of crossings are allowed at the boundary Fi (the 
number of crossings at the other boundary r2 must be 
same modulo 2). Ng = 2^~^ equals the total number of 
valid boundary conditions {qi} G a in each parity sector. 

We calculate entanglement entropy using the reduced 



density matrix, 
latcd, 



Here, 



trB|^)(4'| is readily calcu- 
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1 



2N„ ^ 

{qi}& even 
,A 



-2Rcal(c5oCoi) 
1 

'2n: 



[(IcooP 

.1) 



Coi 



I* 



{«i},Ol 



ground state, the electric charges need to be annihilated 
with each other. One way to do this is to wrap the open 
string O parallel to x around the cycle of the torus. O 
becomes a closed loop C, yet this changes the parity of 
electric field winding number along x. We define the elec- 
tric charge loop operator that insert an additional electric 
.fittld in the x{y) direction by the above procedure as a Z2 
{"Jil^o'sJlpctric flux insertion operator TxiTy). 



{gjeodd 



Clol 



ICllI 



*4},o>(*4},oi 



-2Rcal(c*ioCn) (|* 



{91} 



I* 



i>«},ol 



2Na 



[icoo + coiH* 



A 

{«}.+ 



^ I 



even 



I Coo 



coiPl*^ 



{qi}: 



2N„ ^ L 



cio + ciir|* 



-)(* 



^ I 



{gj6 odd 



ClO ■ 



leob) 

16.) 

ICal) 



(CI) 



There is also a magnetic field, which determines the 
phase of the electric charge as it moves. In particular, 
when there is a magnetic field along the y direction of 
the torus of 1(0) total flux(mo(i2), the electric charge 
picks up a — (+) sign traveling around the loop around 
the X direction, and similarly for the magnetic field along 
the X direction. Denoting the insertion operator of such 
Zi magnetic flux as Fy and , the loop operators of the 
magnetic charge (vison), we have, 



Here 1*^^ , , ) = ^{V^j- ^ ^ ±'^j- ^ M and hold 

the orthogonal condition. 

From the above expression, it immediately follows that 
the Renyi entanglement entropy Sn is given by EgnfTBl 




Nq YPPiT 



= 1 



lo 



1—71 ■' 

3 = 1 

= L\og2- ^21og2 + -l^log^p;^j 
where pj are defined in Eqn. 1191 

To understand the nature of the corresponding MES in 
Eqn. 1201 we first discuss the quasi-particlc excitations of 
the Toric code model. Imagine acting a string operator 
defined on the links of the lattice 



Now W^'{0)\vekCx) is an excited states and still an 
eigenstate of As and Bp , with As = — 1 at the two ends of 
O. We may regard them as electric charge quasi-particles 
that cost a finite energy to create and the string con- 
necting them as an electric field line. To return to the 



TxFy 



Fx Fx Ty 

They suggest that Tx{Ty) is the magnetic fiux measur- 
ing operator in the y{x) direction and Fx{Fy) is the elec- 
tric flux measuring operator in the y{x) direction. Note 
that both electric and magnetic flux are defined mod- 
ulo 2 in correspondence with the Z2 gauge theory. After 
simple algebra. 



Fv\U) 
Fx\U) 



(-l)IU) 
(-I)'IU) 



(C2) 



Compare Eqn. ICll and I C2I with Eqn. [20l we arrive at 
the conclusions listed in Table HI 



Appendix D: Modular Transformations 

The S and U matrices describe the action of modular 
transformations on the degenerate ground states of the 
topological quantum field theory on a torus. For Abelian 
phases, the ij'th entry of the S matrix corresponds to 
the phase the i'th quasi-particle acquires when it encir- 
cles the j'th quasi-particle. The lA matrix is diagonal and 
the li'th entry corresponds to the phase the i'th quasi- 
particle acquires when it is exchanged with an identical 
one. Let us first review the geometric meaning of these 
transformations. Labeling our system by complex coor- 
dinates z = x + iy, the torus may be defined by the 
periodicity of wi and lu2 along the two directions e"i and 
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62 (need not to be orthogonal) i.e. 
Now consider a transformation 



Z = Z + iUl 



UJ2 



a b 
c d 



UJ2 



L02- 



(Dl) 



where a, b,c,d(z TL. Since our system hves on a lattice, 
the inverse of the above matrix should again have integer 
components, hence the determinant ad— be = 1. One can 
show that matrices with these properties form a group, 
called 5*^(2, Z). Interestingly, all the elements in this 
group can be obtained by a successive application of the 
following two generators of S'L(2,Z): 



S = 



This transformation corresponds 



1 

-1 

to uJi — > UJ2 and UJ2 — > — cji and therefore, for a 
square geometry corresponds to rotation of the sys- 
tem by 90°. 



U = 



Under this transformation uji 



1 1 
1 ^ 

uj[ = cji + UJ2 and 0^2 — > = i^2- Consider a 
loop on the torus with winding numbers ni and 
712 along uji and uj2 directions. By definition of 
the U transformation, the winding numbers in the 
transformed basis: 



niCJi + 712012 

ni(uj[ — 0^2) + 'T-2'^2 



This matrix belongs to the group SL{2,Z) and simple 
algebra shows that i?27r/3 = US. One may also check 
that -R27r/3 = 1 as one might expect. Therefore, knowing 



the action of R 
matrix. 



27I-/3 



on the MESs would lead to the US 



Appendix E: Modular matrices of Z2 gauge theory 
by transforming minimum entropy states 



Let's study the action of modular transformation on 
the MESs ISq) for the Z2 gauge theory in Sec. Ill CI and 
compare the resulting modular matrices with the known 
results. 

First consider a 7r/2 rotation symmetric square sample. 
Under tt/2 rotation, \^ab) — >■ \£,ba)- According to Eqnl20l 
the transformation for the MESs for cuts along y: 



^1)^ 2(|Si) 



^(|Sr>- 



^4>-> 2(1^1) 



IS2) + IS3) + IS4)) 
IS2) - IS3) - |S4» 
S2) + IS3) - |S4» 
S2> - IS3) + |S4» 



where 



ni and 773 =712 — 711 are the winding 



numbers along the and W2 directions. 

The transformation properties of the resulting MESs 
under modular transformations would yield the desired 
S and U matrices. Further, for a symmetry transfor- 



mation of F{S, U) on (cji, 



t^2j 



the corresponding mod- 



ular transformation on MESs would yield the modular 
F{S,U) matrix. 

In the main text, we have obtained S and hi matrices 
for the toric code model from the action of these transfor- 
mations on the basis states |^ab). Wc now show that one 
can also obtain the IAS matrix by studying the action of 
27r/3 rotation R2n/3 on the MESs (provided that i?2ir/3 is 
symmetry of the model). To see this, consider a triangu- 
lar lattice that is defined by two lattice vectors (complex 
numbers) wi,a;2 with oji = (1,0) and L1J2 = (1/2, •\/3/2). 
The transformation of our interest is the transformation 
of LUi,LU2 under R2tt/3 rotation: wi — >■ = — cji +i02 and 
a;2 — > ^2 = —LOi. Therefore, one can write the R2tt/3~ 
matrix 



R 



27r/3 



-1 1 

-1 



(D2) 



Hence, the modular S matrix is given by 




S 



This is exactly what one expects from the topolog- 
ical quantum field theory corresponding to the zero 
correlation length deconfined-confined Z2 gauge theory. 
There are four flavors of quasi-particles in the spectrum: 
l,m,e,em, as we have shown in Table HI The electric 
charge e and magnetic charge (vison) m both have self- 
statistics of a boson and pick up a phase of tt when they 
encircle each other (and as a corollary, the same phase 
when they encircle em). By studying 5, one gets the self 
and mutual statistics for quasi-particles encircling each 
other. 



In Sec IIIIBI we further show that symmetry is not re- 
quired to determine the S matrix. In Eqn. [20] we have 
shown the MESs for cuts along W2 = y direction: 
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(l^oo) 

(leoo) 



IS3) = ^(l^io) 



F4 



^(lao) 



leoi)) 

leoi)) 

ICii)) 

lai)) 



where are undetermined phases for MESs jS^). The 
unitary matrix Ui connecting the W2 MESs and the elec- 
tric flux states: 



Ui 



V2 



V 



(El) 



On the other hand, it is straightforward to verify that 
for loops along w'2 — —x + y direction, which satisfies our 
requirement Ean l231 the corresponding MESs: 



This leads to the following modular S matrix: 




which is indeed the correct result for Z2 toric code. 

Now consider the transformation corresponding to U 
matrix as described in Appendix |D1 where n'^ = ni and 
71-2 — 712 — ni are the winding numbers along the uj[ 
and UJ2 directions. Using this expression and Eanl20l the 
transformation for MESs from W2 cut to W2 cut: 



IS2) ^ IS2) 

IS4) ^ -IS4) 





= ^(leoo) 




S'2> 


piv'2 


- iCii)) 


%) 




+ ICio)) 




piv'i 


- ICio)) 



This leads to the following modular U matrix: 



again Lp\ are undetermined phases for MESs The 
unitary matrix U2 connecting the MESs and the elec- 
tric flux states: 



V 



(E2) 



Combining Eqn. lEll and IE2[ we can write down the 
modular S matrix as: 



S ^ U2^Ui 

/ g'i(Vl-Vi) gS'(V2-ip'i) g'i(V3-vi) _g«(V4-vi) \ 

1 g^CVl-vSi) gj(V2-iP2) _g«('l=3-'i=2) g*('l=4-V2) 

~ 2 g^lVl-Vs) _g«(V2-¥'3) g^CVS-Vs) gi(V4-V3) 

y gi(vi-v4) —Qi(v2-v'i) _gi(v3-v4) _gj(v4-v4) y 

To ensure the existence of an identity particle in accord 
with the first row and column, we impose the conditions: 



U 



10 
10 
10 
-1 



Again, this is what is expected from the Z2 gauge the- 
ory. The sign of — 1 on the last entry of the diagonal cor- 
responds to the fermionic self statistics of the em while 
the positive signs correspond to the bosonic self statistics 
of 1, e and m particles. 

To sec a more generic example to derive the lA matrix 
from rotation symmetry, we first define the toric code on 
a triangular lattice, with system dimensions such that 
the 27r/3 rotation is a symmetry of the system. The 
Hamiltonian is same as Eq. [THwith the star 's' denoting 
six links emanating from a vertex while the plaquctte 
'p' now involves three links. We again denote the four 
degenerate ground states on a torus as \S,ab) with a,b = 
0, 1 denoting the parity of electric field along the non- 
contractible cycles. The relation between the MESs 
and the states \^ab) remains unchanged (Eqn. [20)) . The 
calculation for the transformation under 27r/3 proceeds 
analogously to that for tt/2 rotation and one finds: 



^27r/3l^Oo) 
-R27r/3l^0l) 
-R27r/3l^lo) 
^27r/3l6l) 



ICoo) 
ICio) 

l^ll) 

leoi) 



17 



Translating the action of i?27r/3 on the states |^ct) to tains 
that on states IS™), one finds 



/fill 




U = 



10 
10 
10 
-1 



Combining the expression and the S matrix, one ob- 



as expected. 
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